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Abstract 

We show that topological properties of minimal Dirac sheets as well as of 
currents lines characterize the phases unambiguously. We obtain the minimal 
sheets reliably by a suitable simulated-annealing procedure. 



1. Introduction 

We investigate compact U(l) lattice gauge theory in 4 dimensions with the action 
[1] _ 

s = /3 E (l-cose^ + A^lM^I 

{1>V,X p,x 

where M P:X = t pailv {Q ^ x+a - A1 ^)/47r. The physical flux G [-7r,7r) is given 

by [2] Q^j; = Q^j; + 2Trn fll/ ,x- We consider periodic boundary conditions (apart 
from a final remark on our results for open boundary conditions). 

The strength of the first order transition decreases with A, the transition ultimately 
getting of second order [3]. We use this to set up a very efficient algorithm [3, 4] for 
the Monte Carlo simulations in which A becomes a dynamical variable. 
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We argue that analyzing topological structures of the configurations (rather than 
only measuring global observables) is a promising tool. We apply it to monopole 
currents and Dirac plaquettes. 

For loops the topological characterization is straightforward. For the actually 
occurring networks of monopole currents we have recently been able to present 
a mathematically sound characterization [3]. This is not only necessary for the 
unambiguous identification of physical features but also for the computer analysis 
of huge networks. 

Here the analysis of configurations is extended to Dirac sheets and the fact that 
appropriately specified topological structures signal the phases is confirmed in more 
detail. In conclusion the general principle, not relying on the particular boundary 
conditions, is pointed out. 

2. Dual-lattice structures 

On the dual lattice for the current J PjX = M PjX+p related to links the conservation 
law J2 P (Jp,x — Jp,x- P ) = holds. The definition of current lines is such that for 
Jp,x = there is no line on the link, for J PjX = ±1 there is one line, and for J PjX = ±2 
there are two lines, in positive or negative direction, respectively. Networks of 
currents are connected sets of current lines. For a network N disconnected from 
the rest one gets the net current flow / with components = ^2 X x x Jfi 3 ,x for 

The Dirac string content of the plaquettes on the dual lattice is described by 
P P a,x = -^pa^n^x+p+a for which the field equation E^fiW ~ P P a,x-a) = J P ,x 
holds. Dirac plaquettes then are defined such that for p pa , x = there is none, for 
p pa ,x = ±1 there is one, and for p pa , x = i2 there are two of them. 

Dirac sheets are formed by connecting Dirac plaquettes with common edge and 
appropriate orientation so that J PjX = except at the boundaries of the sheets. 
Dirac sheet structures are not gauge invariant. However, they belong to equivalence 
classes which cannot be deformed into each other by gauge transformations. For 
given boundaries there are topologically distinct possibilities for the sheets. Thus 
the equivalence classes of sheets carry more information than the current networks. 

To represent the equivalence classes of sheets by their members with minimal area 
we minimize the number of Dirac plaquettes by gauge transformations. We find that 
careful annealing is necessary to get the minimal structures reliably. For the anne- 
aling we use Metropolis sweeps for P(tj) ~ exp(— arj), where rj = Y^x,p,>u l^^l, with 
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random local gauge transformations keeping the angles in [— 7r,7r). In subsequent 
simulations we increase the parameter a in appropriate steps. 

3. Topological analysis 

The elements of the fundamental homotopy group 7i"i(X, b) of a space X with 
base point b are equivalence classes of paths starting and ending at b which can be 
deformed continuously into each other. Its generators may be obtained embedding a 
sufficiently dense network N into X and performing suitable transformations which 
preserve homotopy. We use the observation that by such a method only the gene- 
rators of a subgroup are obtained if a given network N does not wrap around in all 
directions. This provides an unambiguous characterization of networks. 

Choosing one vertex point of N to be the base point b and considering all paths 
which start and end at b we note that a mapping which shrinks one edge to zero 
length preserves the homotopy of all of these paths. Therefore, by a sequence of 
such mappings we can shift all other vertices to b to get a bouquet of paths starting 
and ending at b without changing the group content. 

Describing a path by a vector which is the sum of oriented steps along the path, 
for N with K vertices and K\ edges the bouquet has K = K\ — K + 1 loops 
represented by vectors Si with i = 1, . . . , K and j-th component Sij = WijLj where 
Lj is the lattice size. The bouquet matrix Wij then is to be analyzed with respect 
to the content of generators of 7Ti(T 4 , b) = Z 4 . 

Current networks have the additional properties of path orientation and current 
conservation. The maps reducing the bouquet matrix w^, therefore, in addition 
to homotopy have to respect current conservation. This leads to a modified Gauss 
elimination procedure within which adding of a row to another one requires to 
subtract it simultaneously from a further row. In this way one arrives at the minimal 
form with rows ai, . . ., a r , t, 0,. . ., where r < 4. The significance of this gets obvious 
switching to the pair form with rows ai, — ai, . . ., a r , — a r , /, 0, . . ., which exhibits 
the relation to the net current flow / explicitely. It is to be noted that / ^ occurs 
only in very rare cases [3] . The number of independent pairs determines the number 
of nontrivial directions. 

For the topological analysis of minimal Dirac sheets we use the network of plaquet- 
tes as a sufficiently dense auxiliary network. The bouquet matrix then is obtained 
as described before. The reduction of the bouquet matrix is simpler here because 
usual Gauss elimination applies, which gives the minimal form with rows ai, . . ., 
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Figure 1: Probability -P ne t for hot and cold phase as function of A for lattice sizes 
8 4 (circles) and!6 4 (crosses). 



a, n 0,. . ., 0. Now the number of independent vectors corresponds to the number of 
nontrivial directions. 



4. Results 

Figure 1 shows our results for the probability P net to find a current network which 
is nontrivial in four directions as function of A in the transition region. For larger L 
or negative A, where the peaks of the energy distribution related to the phases get 
well separated [3, 4], P ne i is seen to be 1 for the hot (confining) phase as compared 
to for the cold (Coulomb) phase. We thus have the unambiguous topological 
characterization by the existence of a nontrivial network in the hot phase and its 
absence in the cold phase. 

Our analysis of minimal Dirac sheets leads to analogous results as that of current 
networks. The hot phase is characterized by the existence of a topologically non- 
trivial sheet and the cold phase by its absence. Because, given a nontrivial current 
network, the sheets could be trivial or nontrivial, this answers the question about 
the location of the sheets. 

We construct minimal Dirac sheets by making first the connections where only two 
Dirac plaquettes meet at an edge with appropriate orientation. Then we consider 
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Figure 2: Size of largest current network (circles, links) and largest minimal Dirac 
sheet (crosses, plaquettes) in hot and cold phase as function of A on 8 4 lattice. 



the places where more than two plaquettes meet at an edge in such a way that locally 
there is more than one possibility to make a connection. We find that connecting all 
and connecting none of the plaquettes at those places makes very little difference in 
the results. Thus we get the further remarkable observation that the minimal sheets 
do not intersect much. 

In Figure 2 we present the sizes of the largest current network and of the largest 
minimal Dirac sheet. Obviously the sheets give a more sensitive signal for the phases. 
In Figure 3 we demonstrate that our probability P net is an order parameter superior 
to n max /n tot , the relative size of the largest network, advocated in [5]. 

In order to consider our results from a more general point of view we note that 
a topologically nontrivial structure on a finite lattice with periodic boundary con- 
ditions corresponds to an infinite structure on an infinite lattice. This suggests to 
characterize the hot phase more generally by the existence of an infinite network (or 
sheet) and the cold phase by its absence, where on finite lattices "infinite" is to be 
interpreted according to the particular boundary conditions. To test this characte- 
rization we have also performed simulations with open boundary conditions. They 
show that outside the transition region (which now is extended due to finite size 
effects) there is still a clear topological signal for the phases provided that in the 
analysis the prescription "nontrivial in all directions" is replaced by "touching the 
boundaries in all directions". This supports the indicated percolation-type picture. 
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Figure 3: Order parameters P net and n max /n tot for A = as function of f3 on 8 4 
lattice. 
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